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Abstract 


We investigate the stability of an inhomogeneous chiral condensed phase against low energy 
fluctuations about a spatially modulated order parameter. This phase corresponds to the so-called 
dual chiral density wave in the context of quark matter, where the chiral condensate is spatially 
modulated with a finite wavevector in a single direction. Prom the symmetry viewpoint, the phase 
realizes a locking of flavor and translational symmetries. Starting with a Landau-Ginzburg-Wilson 
effective Lagrangian, we find that the associated Nambu-Goldstone modes, whose dispersion rela¬ 
tions are spatially anisotropic and soft in the direction normal to the wavevector of the modulation, 
wash out the long-range order at finite temperatures, but support algebraically decaying long-range 
correlations. This implies that the phase can exhibit a quasi-one-dimensional order as in liquid 
crystals. 
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I. INTRODUCTION 


Unveiling the phase diagram of quantum chromodynamics (QCD) is among the most fun¬ 
damental issues in quark-hadron physics. So far, considerable theoretical and experimental 
efforts have been devoted to exploring the QCD phase diagram (for reviews see Ref. [l|); 
the properties of the high temperature regime are studied experimentally in ultra-relativistic 
heavy-ion collisions, and in ab-initio lattice QCD simulations. The latter are subject to tech¬ 
nical difficulties at non-zero net baryon-number density, the so-called sign problem. In the 
near future, data at lower beam energies, relevant for the exploration of the phase diagram 
at non-vanishing baryon density, will be forthcoming. In order to exploit this opportunity in 
an optimal way, it is necessary to hnd appropriate observables for deciphering the properties 
of dense and moderately hot matter in such collisions [^. 

In recent theoretical studies of QCD at hnite temperature and density, various inhomoge¬ 


neous chiral condensed phases have been proposed (for a recent review see Ref. [3|). These 
studies suggest that the conventional QCD phase diagram should be redrawn. Indeed, it is 
possible that the phase structure at high net baryon densities and moderate temperatures is 
modihed considerably by the presence of inhomogeneous phases. Thus, the region of the chi¬ 
ral transition may be extended and the order of the phase transition may change. These fea¬ 
tures are gleaned primarily from mean-held calculations in the Nambu-Jona-Lasinio (NJL) 
and quark-meson (QM) models ^la], and the Dyson-Schwinger approach to dense QCI)_j6[. 
It also interesting to note that within the Gor’kov approach to chiral effective models [7,|Jit 
is found that the QCD critical endpoint is a Lifshitz point, where the normal, homogeneous, 
and inhomogeneous chiral condensed phases meet. In the large Nc approach to dense QCD, 
early studies suggested the emergence of the so-called chiral density wave {q, 10|, while in the 
context of quarkyonic matter ll| another inhomogeneous phase, the so-called quarkyonic 


chiral spiral, was discussed 


12 |. 


The inhomogeneous chiral condensed phases mentioned above correspond to a one¬ 
dimensional modulation embedded in three spatial dimensions. Some of these structures 
are based on extrapolations from analytic solutions obtained for purely (1-1-1)-dimensional 
systems j^. Possible extensions to higher-dimensional modulations have been studied, with 
the result that the one-dimensional modulation tends to be favored close to the Lifshitz 
point and/or at zero temperature jl^ . 
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Let us start by classifying the modulations for inhomogeneous phases, according to the 
convention employed in condensed matter physics. There are basically two types of one¬ 
dimensional modulations: one is of the Fulde-Ferrell (FF) type [^, characterized by modu¬ 
lations of the phase of a complex order parameter with constant amplitude, while the other 
is of the Larkin-Ovchinnikov (LO) type [^, where by contrast only the amplitude is mod¬ 
ulated. The FF type includes the dual chiral density wave (DCDW) ^ and the qu arky onic 


chiral spiral 


12|. On the other hand, the chiral density wave (a plane wave) jol, 10| and 


periodic domain walls [5| are of the LO type. In the present paper we focus on the DCDW, 
which is of FF type. The DCDW is characterized by modulated scalar and pseudoscalar 
condensates with a constant amplitude A and a wavenumber g jd]. 


= A cos qz, {4’i'y5T3'ip) = Asingz, 


( 1 ) 


where '0 is the quark held for two havors, and Ts a Pauli matrix diagonal in the isospin 
space. This conhguration is akin to avr® condensation, obtained in neutron matter within 


18| . and is thus expected to smoothly connect between nuclear and quark 


the sigma model 
matter. 

Most studies of inhomogeneous chiral condensed phases so far are based on mean-held 
calculations. Thus, ehects of thermal and quantum huctuations have yet to be studied. In 
the context of pion condensation, the stability of modulated condensates against thermal 


huctuations has been studied on the basis of Landau-Peierls arguments 


19| . It was found 


that in such systems there is no true long-range order with a non-vanishing order param¬ 
eter [^. Instead, such systems can develop a quasi-ordered one-dimensional condensate, 
with correlation functions that decay algebraically in space.* In this paper, we investigate 
the stability of the DCDW phase against low-energy huctuations of Nambu-Goldsto ne ( NG) 
modes associated with the spontaneous symmetry breaking, along the lines of Ref. 2C |. 

The paper is organized as follows. In the next section, we construct a (3-1-1)-dimensional 
Landau-Ginzburg-Wilson ehective Lagrangian for general order parameters of the chiral 
condensate, which are allowed to be spacetime dependent, and then apply the formalism to 
the DGDW phase. In Sec. IIIIl we discuss the symmetry breaking pattern and the corre¬ 
sponding NG modes in the DGDW phase. We also present the dispersion relations for these 


* A similar state is found, e.g., in smectic liquid crystals 2l|, [22|. See also Refs. 23|,l2^ for a corresponding 
discussion on FFLO superconductors/superfluids. 
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low-energy collective excitation modes, by introducing fluctuations such as amplitudons and 
phonons/phasons (NG modes) on the ground state of the DCDW. In Sec. HVl we investigate 
the impact of low-energy fluctuations on the order parameter by evaluating the long-range 
correlation functions. Finally, section |V] is devoted to summary and outlook. 


II. LANDAU-GINZBURG-WILSON EFFECTIVE LAGRANGIAN 


We start by introducing a 2 x2 matrix held M as the (1/2,1/2) representation of the chiral 
SU{2)l X SU{2)r symmetry. In the following, we shall use the fact that SU{2)l x SU{2)ji 
is isomorphic to the four-dimensional rotation group 0(4). The matrix M can be expressed 
in terms of the right and left handed quark fields 25|, Mij = Under the 

transformations Ur^l = exp —i ± • r/2 , where r is the isospin Pauli matrix and a 
and (3 are two three-dimensional vector parameters, the matrix M transforms (to leading 
order in a and /3) as M^j Mij -\- • [r, M]^ - + ^j3 ■ {r, M}^y With the parameterization 

M = a + m -f, one finds the corresponding transformation laws for a and if: a —)■ a — /9 • if 
and if—^-if — axif-|-/3 cr. Thus, the rotation with a corresponds to the vector (isospin) 
rotation while that with /3 to an axial vector (chiral or axial isospin) rotation, respectively. 
We can then introduce a four-component composite held (fP" = (cr, if), which transforms as 
a four-dimensional vector under 0(4) rotations. 

Now, we construct an effective Lagrangian density C with 0(4) symmetry in terms of 0 
and its derivatives,^ 


c = C29o0-9o0-V, (2) 

V = 020 • 0 -f 04^1 (0 • 0)^ -f a4^2V0 • V0 -f a6,lV^0 • V^0 

+06,2 (V0 • V0) (0 • 0) + 06,3 (0'0)^ + 06,4 (0 ' V0)^ , (3) 


Note that since the Lagrangian describes the excitations of a medium, we do not assume 
Lorentz invariance. The potential term V is expanded up to sixth order in powers of the 
held and fourth order in its derivatives, as required for stability of the inhomogeneous phase 
at the mean-held level. 


'1^ For non-zero isospin charges, a term of the form ^eijkiRijQij is added to the Lagrangian. Here €ijki 
(ei 234 = 1) is the antisymmetric tensor, and Rij = —Rji the chemical potential for the isospin density 
Qij — d0ido0_7 0jdo0i). 
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Hereafter we set C 2 = 1 for simplicity. The expansion coefficients a* j in V can be evalnated 


within effective chiral models, like the NJL 


Q and QM js, 8 


8 | models. In the former, one 


finds the following relations among them 04^1 = 04^2 and (og,!, 00 , 2 , ^ 6 , 4 ) = (1/2, 3, 2) 00 , 3 . 


III. LOW ENERGY EFFECTIVE MODES IN DCDW PHASE 

We consider an inhomogeneons time independent chiral condensate of the DCDW type, 

0Q = A (cos gz, 0, 0, singz), (4) 

where A is a constant amplitnde corresponding to and q is the wavennmber 

of modnlation in the 2 ; direction. The valnes of A and q are determined by minimizing the 
potential term of the Lagrangian. For the condensate (po, the potential term reads 

V(0o) = 0 ^ 2 ^^ ® 4 ,iA^ + a 4 ^ 2 g^A^ + a 6 ,ig"^A^ + oggg^A^ + ag^ 3 A^. (5) 

Stability of the inhomogeneons phase is gnaranteed by 

Og^ > 0 (or 06^3 > 0), and ag^iOg^s — ag2/4 > 0. (6) 

The stationary conditions for q and A, = 0, yield 

2gA^ ( 04,2 + 2ag^ig^ + ag^2A^) = 0, (7) 

2A 1^02 + 04 , 2 *?^ + + 2 ( 04^1 + ciggg^) A^ + 3ag^3A^j = 0. ( 8 ) 

They admit three types of solntions: 

(a) Normal phase: g = A = 0, 

(b) Homogeneons chiral condensed phase: g = 0, A 7 ^ 0, 

(c) Inhomogeneons chiral condensed phase: g^ = — ( 04,2 + ag, 2 A^) / 2 ag 4 , A 7 ^ 0. 

For a given set of coefficients, the phase with the lowest energy, is realized on the clas¬ 
sical level. The coefficients implicitly depend on the medinm, and are thns fnnctions of 
thermodynamic variables, like temperatnre and chemical potentials. 
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A. Symmetry breaking and Nambu-Goldstone modes in DCDW phase 


In the DCDW phase with non-vanishing A and q, the SU{2)l x SU{2)r chiral symmetry, 
as well as the translational invariance in the 2 ; direction and the symmetry under rotations 
about the x and y axes are spontaneously broken. To see the symmetry breaking pattern 
explicitly, we hrst perform inhnitesimal transformations corresponding to a spatial transla¬ 
tion in the 2 ; direction with a displacement parameter s and a chiral rotation through the 
angles a and /3: 


<^0 00 + ^ 


^ — {sq /ds) sin qz ^ 
01 cos qz — 0^2 sin qz 
02 cos qz + ai sin qz 
V {sq + 03) cos qz ) 


(9) 


The form of the hrst and the fourth components implies that 0o is invariant under a simul¬ 
taneous spatial translation and axial isospin rotation about the z axis through the angle 03 , 
if gs - 1 - 03 = 0. Thus, in the DCDW phase a locking of axial isospin rotations with transla¬ 
tions is realized. In terms of symmetry generators, there are two unique orthogonal linear 
combinations of s and 03 ; one corresponding to a broken generator, the other to an unbro¬ 
ken one. Consequently, the corresponding NG mode is generated by a transformation with 
gs-|-03 7 ^ 0, i.e., by one whose generator is broken in the DCDW phase. In the following, we 
use 03 = 03(f,a;) and s = 0 to generate the NG mode associated with the broken generator. 
Similar arguments for the spontaneous breakdown of internal and spacetime symmetries are 


given in Refs. 


26 |, 


27|. 


The rotations through 0i and Q !2 generate variations in the second component in ll9ll . 
However, the corresponding NG modes are linearly dependent in the sense discussed in 


Ref. 


28l |. In case of a vanishing wavenumber g = 0, only 0i is relevant. Thus, we generate 


the corresponding NG mode using 0i = 0i(t,x) and 02 = 0. Analogous arguments can be 
applied the third component in ([9|), thus eliminating ai in favour of 02 . 

Spatial rotations about the a:-axis by an angle 6^ yields a transformation, which is nonuni¬ 
form in space: 2 ; —)■ ^^cos^a, -|- ysinO^. Similarly, the rotations about the y-axis by 9y yields 
the analogous transformation. However, the corresponding NG modes and those generated 
by translations are also linearly dependent 


We conclude that, although there are eight broken generators for internal and space time 
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symmetries in the DCDW phase, only three independent NG modes remain. These can be 
chosen as the axial isospin rotations generated hy (3 = (3{t,x). 


B. Low energy collective excitations 


= (A + h) 


= iA + 5)U{f3i 


( cos {qz) ^ 
0 
0 

y sin (qz) J 


( 10 ) 


We now consider a general fluctuation in the DCDW phase: 

^ cos {qz + /Sa) cos (32 cos /9i ^ 
cos {qz + /Ss) cos (32 sin (3i 
cos {qz + /ds) sin/92 
\ sin {qz + /^s) J 

Here 6 is the amplitude fluctuation, the parameters (3 = {(3i, (32, (3s} specifles a rota¬ 
tion in the four-dimensional space spanned by the a and vf fields. Finally, U{f3i=i^2,3) '■= 
gi/3iLigi^2r2g*/33i'3^ where Ti^ 2,3 are the 0(4) (axial isospin) generators 251]. This parametriza- 
tion clearly shows that the displacement in the ^ direction is equivalent to a chiral rotation 
through /da. To leading order in the fluctuations, (ITOD yields 

^ —/9a sin qz \ 

(3i cos qz 
(32 cos qz 
\ (3s cos qz J 

which exhibits the fluctuation of the amplitude in addition to the fluctuations corresponding 
to the NG modes. In the following we consider local fluctuations, promoting the parameters 
5 and (3 to fields 5{x) and (3{x), where we use the compact notation x = {t,x}. 

Plugging the above parametrization into the Lagrangian, we can systematically derive a 
low energy effective theory by expanding in powers of the fluctuation fields 6 and (3. Up to 
the second order in the fields, the Lagrangian L = j d^x£ reads 


4> = {1 +6)(j)o + A 


+ 0(l3lSl3uS^), 


( 11 ) 


where 


C = {doSf + A\do(3uf + A\do(3sf - (V^ -F + Vg ), 

= M252 + a6,4A2(Vh)2 + 4a6,ig2(V,5)2 + a6,i(V2h)2, 

V6/3 = 4gA [a6,2A25 - 2a6,iV2h] V,/3s, 

V/3 = a6,iA2 + aepA^ [(v'/ds)' + 4q^ (V./ds)' 


( 12 ) 

(13) 

(14) 

(15) 
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with the mass term = 4 ( 04 ^+ 06 , 2 ?^)^^ + 12 a 6 , 3 A^ and the transverse held (3jj = (3t cos qz 
where (3t = {/9i, /52}- In the above eqnations the stationary condition a 4 , 2 +n 6 , 2 A^+ 2 g^a 6 ,i = 
0 has been nsed. For details of the derivation, see Appendix El 

In order to investigate the thermodynamics of the system, we now move to Enclidean 
space: t —)• —it with the period 0 < r < /9 where jS = 1/T is the inverse temperatnre. For 
Ganssian hnctnations, we obtain the Enclidean action in Fonrier space (for details we refer 
the reader to Appendix |B|) . 



where we have nsed the shorthand notation; ^ dk = T'^^J and k = {ujn,k) with 
the Matsnbara freqnency Un = 27mT(= iu). The inverse propagators in the above matrix 
notation are given by 


g{k) = 2iq oggA^ + 2a6,iP k^, 
- 06,1 dq^kl + (k‘^Y 


and G{k) = u — og,! ( k + 2qk, 


(17) 


We note that for a non-vanishing wavennmber q, the 6 and /dg hnctnations mix. Moreover, 
transverse hnctnations Pt with diherent momenta k and k + 2qz mix, owing to the scattering 
of hnctnations oh the backgronnd modnlation. 

The determinant of the hrst matrix, S'^^(fc)S'(C^(/c) -|- g^{k) = 0, yields the dispersion 
relations of the normal modes involving 5 and /ds, 

2 ' 


u:_y. ~ M -\- ag,i 


u^+kl + ( k‘ 


-I- k + Ak kl + Bkl, 


uj_ ~ ag,i 


ul_kl + (k 


- Ak^ki - Bkt, 


where wG = 


- 4^^ (1 ± = 4 q^A^aQ ,2 - A^afj,2(^A) /m, and 


(18) 
(19) 
B = 


— (2gA^ag.2) /M®. Note that in the massless mode a;_, the dependence on the transverse 
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momentum is subleading, 0{k^)} Consequently, the transverse fluctuations are softer that 
the longitudinal ones. 


Similarly, equating the determinant of the second matrix to zero, S'q ^{k)SQ ^{k) — G‘^{k) 
0 , we obtain the dispersion relation for Pt, 





( 20 ) 


The second term with the negative sign is a higher order correction of 0{k^), stemming from 
interactions with the background modulation. This term is irrelevant for the effects of low 
energy fluctuations and is therefore dropped in the following discussion. 

IV. IMPACTS OF LOW ENERGY FLUCTUATIONS 

At low temperatures, the low energy fluctuations about the calssical DCDW state domi¬ 
nate. We evaluate the contribution of Gaussian fluctuations to the partition function 



( 21 ) 


Higher-order derivative corrections are dropped, with the assumption that fluctuations at 
energies above some cutoff A have been integrated out in the effective Lagrangian ([2]), which 


then involves only the low-energy fluctuations, 6 and /3, explicitly. 

We first explore the impact of low energy fluctuations on the order parameter. 


((A + i)f/(ft)W = A{U0,)^„} + {SU0,)M, 


( 22 ) 


where we use the compact notation {■■■) = / [T>5] VA(3 ■■■e ^^jZ. In the Gaussian 


^ The sign of u: 


2 


'Z — 


in bj- is always positive if 04.1 > 0 , which generally holds if the transition is second order. 
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approximation, the two contributions to the expectation value reduce to 




and 


^ cos(g 2 ;)e ^ 

0 
0 

sin(gz)e“*‘^3^/^ J 


( 23 ) 


(517(A) 0o) - 


^ —sm.{qz){5l3^)e ^ 

0 
0 

\cos{qz){6f33)e~^^3)/‘^ j 

Here the second order fluctuations are given by 

(5(a:)/33(x)) ~ 0, 

.2io 2 f \\ 1 f ^ 

1 a2/o2/ ^\ 1 ^ ^ 

and A^(/d|(x)) ~ 


(24) 


(25) 

(26) 

(27) 


2 J (27r)3 ’ 

where the fluctuations {/3f2 3 {^)) ^-re all logarithmically divergent due to the soft modes in 
the transverse directions. Details of the derivation are given in Appendix O Consequently, 
the low-energy fluctuations wash out the order parameter, i.e., they destroy the off-diagonal 
long-range order. 


{(A + 5)f/(A)A) = 0. 


(28) 


This result implies that a DCDW phase with true long-range order strictly speaking does 
not exist at non-zero temperature. Such a phase may, however, be realized in a modified 
form, with a quasi-long-range order (QLRO), analogous to that in the Berezinsky-Kosterlitz- 
Thouless phase in two-dimensional systems 29| and in smectic liquid crystals 2^. As we 
discuss in the next section, the quasi-long-range order is characterized by a power-law decay 
of the order parameter correlation function. 

At zero temperature, on the other hand, quantum fluctuation are not strong enough to 
break the modulating order. In fact, at T = 0 the second order fluctuations are given by 
the infrared convergent integral A^(/9^2(3)(^)) — 11 obtained by taking the zero 

temperature limit of Eqs. flC6p and flC12p in Appendix 0 
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The results of this section imply that the transition temperature of the true DCDW 
phase is Tdcdw = 0. Now assume that there is a critical temperature Tc > 0, where the 
system becomes unstable with respect to the formation of a state with a modulated order 
parameter. Then a quasi-one-dimensionally ordered phase or a phase with true long-range 
order in two or three dimensions may be realized below T^. [^. To determine which phase 
is preferred, one must, in principle, compare their free energies. Considering the different 
nature of these phases, this is a challenging task. 


A. Long-range correlations 

We now explore the behavior of the correlation functions in the Gaussian approxima¬ 
tion. Since the order parameter is vector-like, we dehne correlation functions among the 
components: 


fij{x) = ((/)i(a;)0*(O)). 

(29) 

These correlation functions are spatially anisotropic owing to the one-dimensional modula¬ 
tion of the background. We compute the dependence of the iso-scalar correlation function 

on The diagonal components which contribute to the scalar channel are of the form: 

fnizz) ~ -^A^cosgze )^>/2^ 

8 

(30) 

where (3^ = (3i{z) — /Sj(0). For details we refer to Appendix iDl 

Similar results are obtained for the other components: 


f22{zz) ~ -^A^ cosgze“^'=i’2,3((/5i )^)/2^ 

8 

(31) 

fssizz) ~ -A^ cos gze“^'=2,3<(/5i )^)/2^ 

(32) 

and fAA{zz) ~ ^A^ cos )^)/2_ 

(33) 


Here the exponents, ((/5j )^), exhibit the following functional form at large z, 

T f d^k 1 — cos kzZ 


((/^M(3))^)/2 - 


2A^J {271] 
T 


3 


In 




2q ’ 


(34) 


where A is an ultraviolet cutoff. 
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Putting it all the together, we obtain the long-range scalar correlation in the z direction, 


{(i){zz) ■ 0*(o)) 


1-2 ( Z 

~ -A cos qz — 

2 \2;o 


-T/To 


(36) 


where Zq = 2g/A^, and Tq = 167ra6^iA^M^_. In the similar manner, we compute the form of 
the long-range correlation function in transverse directions, 

-2T/To 




iAd- 

2 \xo 


(36) 


where xq = A“^, xt is the transverse distance, and the factor 2 in the exponent of xt/xo 
reflects the number of transverse directions. Note that, in contrast to the longitudinal 
direction, there is no modulation of the correlation function in the transverse directions. 

In this section we have shown that quasi-long-range order of the one-dimensional DCDW 
phase feature algebraically decaying correlation functions at large distances. The slow decay 
of the spatial correlations distinguish the quasi-ordered phase from normal or disordered 
phases, characterized by exponential decays. Depending on the experimental resolution and 


hnite size effects, the algebraic correlations can effectively mimic true long-range order 
2l|,l29|. 


20 


V. SUMMARY AND OUTLOOK 


In this paper we have explored the soft modes of an inhomogeneous chiral condensed 
phase with one-dimensional modulation, the DCDW phase. We found that this phase ex¬ 
hibits a flavor-translation locking symmetry and clarihed the counting of Nambu-Goldstone 
modes. The dispersion relations for collective excitations, including the NG modes, were 
derived. The low-energy modes are spatially anisotropic and particularly soft in the direc¬ 
tions transverse to the modulation, owing to the lack of terms quadratic in the transverse 
momentum in the dispersion relations. As in smectic liquid crystals, the absence of such 
terms is a consequence of the symmetry under rotations about any axis orthogonal to the 
modulation direction 


19 


22 |. 


Moreover, we have shown that at non-zero temperatures the DGDW phase exhibits 
a Landau-Peierls instability, i.e., the long range order is destroyed by low-energy (long- 
wavelength) fluctuations of the order parameter. Nevertheless, a phase similar to the smec¬ 
tic phases of liquid crystals, characterized by a quasi-long-range order with algebraically 
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decaying order parameter correlation functions, is possible. Such an “algebraic order” can, 
depending on the conditions, emulate true long-range order. In particular, this would be 
the case, in a hnite systems, where the range of the order-parameter correlations exceeds 
the size of the system. 

The experimental verihcation of “algebraic order” can be challenging. The slow decay of 
the correlations has been observed by light scattering in smectic-A liquid crystals 


neutron scattering in Bragg glass 


30 


211, by 


Berezinsky-Kosterlitz-Thouless type 


and only recently in a two-dimensional system of the 


3ll , by measuring the coherence of photons emitted in 


quasiparticle decays. Whether the quasi-one-dimensionally ordered DCDW phase could be 
observed by an appropriate choice of probes is still an open question. Hence, it would be 
important to systematically explore how the collective modes in the DCDW phase interact 
with external probes such as hadrons (quarks) and photons (gauge helds). 

There are also several theoretical issues, that deserve further study. In particular, it 
is known that inhomogeneous chiral phases are favored in systems with vector-vector type 
interactions, which tend to enhance the size of the inhomogeneous area 32l, and that in 


the presence of an external magnetic held 


temperatures, by topological aspects 


33, 


he FF type phase is stabilized, also at hnite 


3^. Moreover, since two- and t 


iree-dimensional 


condensates with true long-range order are allowed at any temperature 


19| . it would be 


important to compare the free energy of such phases with that of a one-dimensional conden¬ 
sate. It would also be interesting to understand how higher order interactions among the 
collective modes modify the soft modes. These may ahect the Landau-Peierls instability of 
inhomogeneous phases discussed here. 

Finally, the topics discussed here may have an impact on the physics of compact stars. 


It has 
phases 


Deen speculated that various spatially inhomogeneous phases, like nuclear pasta 


35l | and hadron-quark mixed phases 3^, could be realized in the interior of such 


stars. These could have phenomenological implications, allowing, e.g., novel cooling sce¬ 
narios (37|. Thus, it would be interesting to study the properties of inhomogeneous chiral 
condensed phases under conditions relevant for neutron stars in general and compact stars 
with quark cores in particular, i.e., in charge neutral matter in /3 equilibrium but also at 
nonzero isospin density 38| and hnite strangeness 39|. These topics will be considered in 


future works. 
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Appendix A: Effective action of flnctuations 


We derive the effective action of the fluctuations to the second order of held expansion. 
It mostly comes from derivative terms: (V0)^ and (V^0)^. First of all, we rewrite the held 
as 

(t){x) = (A + 6)U {l3i)(t)o{x) = (A + (5)I^(A)0o, (Al) 

where V = US with S = and 0o = (1, 0, 0, 0)"^. 


1. term 

We derive the derivative term (V0)^ up to the second order of huctuation helds: 

V0 • V0 = 0^ + (A + 6)VV-^] [V(5K + (A + 6)VV] 0o 

~ (V(5)2 + 0^6^ + 2A6 (g2 + , (A2) 

where 

0^ 00 = g' + 2gV,03 + + (V03)^ (A3) 

and we have used the stationary condition under which 04 . 2 , Oe.i; and 00.2 terms disappear. 

2. (V^0)^ term 

We next derive the derivative term (V^0)^; 

V20 . V20 = 0^ + AVV-^] + AVV] 00 

= + A(V^(5K-^)VV) 

+A(VV-V2((5l^)) + A2 (vV-ivV), (A4) 

where (■ ■ •) = ■ ■ ■ ^o- We expand each term up to second order of the huctuation helds. 

Hereafter we chop oh terms which will disappear together with 04 . 2 , ag.i, and 05.2 terms 
under the stationary condition, and the constants. Omitting total derivatives, we obtain 

~ {V^S-q^6f + 4q\Vjf, (A5) 

{V"(<5H-')VV) ~ q^d + 4gV2V,<503 - VV,<503, (A6) 

(VV-W) ^ ([(V' + g") {S-^p)] [(V' + g") {pS)]). (A7) 


17 



where ^ 


Appendix B: Propagators 


The following propagators are derived from Eq. m-- 


f {6ik)6*ik)) {6{k)Af3;ik)) \ 

{A^sik)5*ik)) {A^sik)A^m) ) 

^ 1_1_ f So\k) g{k) \ 

2Sso\k)S^\k)+g^{k) y _g^k) S;,\k) ) ’ 

I {A(3,{k)A(3*{k)) {A(3,{k)Al3t{k + 2qz)) \ 

y{Al3i{k + 2qz)Al3*{k)) {Al3i{k + 2qz)Al3*{k + 2qz)) j 

1 1 I So\k + 2qz) -G{k) 

~ 2 S^\k + 2qz)S^\k) - G\k) 1 _G{k) Sk\k) 


(Bl) 


(B2) 


Here note that each component of the latter matrix S{k) corresponds to summing all tree 
diagrams of the absorption/emission processes, for instance, the {1,1} component reads 


= So{k)Y,[G{k)So{k + 2qz)G{k)So{k)f 

n=0 

_ Sk^{k + 2qz) 

~ Sk\k + 2qz)Sk\k)-G^{k) 


(B3) 


Appendix C: Fluctnation effects on the order parameter 


In averaging the order parameter over quadratic fluctuations in Eq. fl2^ . we have used 
following results: 


(-ly 


{cos{qz +(5^)) = cosqz^-^{l3l^) 

n—n ' h 


= COSqze 


n=0 
-(/3|>/2 


(d COS (qz + ks)} = — sin qz 


( /o2n 


-(2n + l)! 


(kt)(2n + l)m 


n=0 


= —sinqz{5(3'j)e 


(Cl) 


(C2) 


and similarly for the other contributions. 
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Also, the following finds useful in evaluating integrals like f d^kT/u^, 



dk^ 


1 

ak^ + b + ck^ 


TT / 1 



(C3) 


where (a, b) = (a, b)/c, = ^ — b > 0, and = f ±A. In general case where a = ci-\-C2k1, 

b = cAti and c = 1, for small kt, the above integral results in as follows: 

TT / 1 1 \ TT / 1 1 

2Ac\B_ B+) ciVVca^f ^/ci 

where 



A 


^ + ^k? 


-kf + 0(kf), and 5|~ 


Cl + C2k^ 


^kf 

Cl t 


(C5) 


Here, in a case that a = Aq^ + b = kf, and c = 1, which leads to J dkzOJZ^, the integral 
can be evaluated as n/2qk‘f, where kt is the momentum in the x-y plane. 

Hereafter we consider the expectation values of second order fluctuations in Eq. fl22l) . 


1. Second order fluctuations for 5 and ds 


We evaluate the second order fluctuations by considering infrared (IR) singularities, and 
in high-temperature and low-energy expansion. 

First of all, the second order fluctuations for (3^ result in as follows: 


c-l 

*^<50 


{k) 


S,{k) ^ 2{A(3,{k)Am)) = ^ 

S6o{.k)SQ [k)+g^[k) 




Ull 


for small \k\, 


^1-^5 1 


+ 


a;| — u'i 



SsA, k) 


^ 3 ( 0 , fc) 


= T^e-^^-^SAk) 

n 


— ul n{u:+)e^+'^ -f (n(a;+) -|- 1) e 
oj\ — uAl 20 ;+ 




bJl — u?_ + 1) e 

^ 0 ;+ — o;^ 2o;_ 

T 

— for small \k\ and o;±/T, 


(C6) 

(C7) 


where o;| = -|- aQA^Ak"^ + 4a6.ig^(fc^)^ + a6.i(fc^)^, and n{x) = — l) ^ the Bose 

distribution function. Here the T = 0 limit of fICbll reads 53 ( 0 , k) ~ 


19 













Secondly, the second order fluctuations for S result in as follows: 


So\k) 


ojI-uI 1 


col -Oj'i 1 


Ss{k) _ 2{6{k)S*{k)) J ^2 _ ^2 + ^2 _ ^2 ^2 _ ^2 

1 ^ 
for small \k\, 




Ss{T,k) = Tj2e-'‘^""Ss{k) = ^ 


co\ — co'p n{co+)e^+'^ + (n(a;+) + 1) e ‘^+'’ 


OJl — 0J_ 

n ~ 

cofj — ul n{cjJ)e^-'^ + (n(a;) + 1) e“ 


2a;_i 


/ )2 , ,2 

CJ I — oj_ 


2cj_ 


T 


SsiQ^k) ~ ^ for small \k\ and a;±/T, 


cut 


(C8) 

(C9) 


where a;| = Og,! 4g^A;^ + , and the T = 0 limit of fICSp reads S'3(0, k) ~ 

Finally, the second order fluctuations for 5 and (3^ result in as follows: 

Ss,{k) = 2{S(k)APm) = , 2^M = - ^-^2 ; • 

^so i^)^o i^)3-g^{k) ul_ — ui \u^ — ul_ cj^ — cot/ 

^53(r,fc) = Tj2e-'‘"’^^Ss3{k) 


9{k) 


col — col_ 


n{cj-)e^-'^ + (n(a;_) + 1) e n(a;+)e‘^+'^ + (n(a;+) + 1) e '^+'’ 


5,3(0, fc) ^ 


CJl — cot 


T T 

uol_ ul 


2co- 

2iqA‘^aQ,2kz T 

AP ZK 


2co4 


for small \k\ and co±/T. 


(CIO) 

(Cll) 


From the above results, the second order fluctuations for S and (3^ in real space read 
= I / ^*^3(0; k)cc^lJ and A(5(a;)/?3(a:)) = ^ / ^5,3(0; k) ~ 0. 


2. Second order fluctuations for /3i and /32 


Here for later convenience, we calculate the imaginary time correlations of mean square 
fluctuation of the /5i^2- 

For the /9i 2 with same momentum. 


Sii{k) 

Sn{r,k) 

Su{0,k) 


^ 2{Af3,,2{k)A/3l,{k)) 


SQ^{k + 2qz) 

~ So\k + 2qz)So\k)-G^k) 

4^ + {nM + 1) , 

ZCOk 


T 

~ —^ for small \k\ and cOk/T, 


(C12) 

(C13) 
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where Uk is the dispersion relation for /3i^2 given by Eq. fl20|) . and fi{k) = j 


-2 


1 for \k\ ~ 0 with := (/c^ — Here the T = 0 limit of flC12p reads S'ii(0, k) ~ 

1 

Similarly, for the /5i^2 with different momentnm, 

- 2(Aft,(.)Aft-,(. + 2,1)) ^ ^ 3 Jg},) _ 

^i2(r, k) = tT e-*^""^i2(fc) = 4^ + (n(a;fc) + 1) , (C14) 


*S'i 2 ( 0 , fc) ~ —9 for small Ifcl and (n^/T, 


-2 


~ 


8g^ a;;;, 

where /r(«) = (l - 

T = 0 limit of fICldp . 

From the above results, the real space fluctuation reads 

A2(Ati,2(a;)) = yf^yfdpdke^^P+^'>^A/3,{p)A/3,{k)) 


(CIS) 


for \k\ ~ 0. Here S'i2(0,A;) ~ ^ 


1 


8(j^ 2wk 


for 


= 2 cos2qz d/cRe(A/?*(A; + 2qz)Af3i{k)) 

+2 sin 2gz dfclm{ A/j* (/c + 2qz)Af3i{k)) + ''^^dk{Af3i{—k)Af3i{k)) 


= 2 cos 2qz 


d^k S2i{0;k) f d^k Sn{0;k) 


+ 


(C16) 


(27r)3 2 J (27r)3 2 

where in the last line the hrst term gives zero due to the odd function, and the second 
diverges logarithmically at fc = 0 at hnite temperature. 


Appendix D: Order parameter correlations 


The following integral is useful in evaluation of order parameter correlation functions at 
long range, 

cos{kzz) TT / e~^-^ e~^+^ 


J_^ ""akl + b + ckt ~ 2cA\ B_ B+ )' ^ ’ 

where = ^ —6 > 0, = |±A, and (a, h) = (a, h)/c. In general case where a = Ci + C 2 k‘f, 

b = c^kf, for small kt, the above integral results in as follows: 


TT 


^—B—Z 


2cH V H_ 


B, 


TT / e 


-J^k 


Cl 




(D2) 
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where 


A = ^+ + 0{kf), and ~ 


Cl + C2k'f 


^kf 

Cl t 


(D3) 


With the above results, we can evaluate the following integral, 


(27r) 


dk,, I dktkt 


1 — cos{kzz) 


/•^ dktkt TT 
L (27r)2ci 

^dktl-e 


(ci + C2k^)kl + c^kf + ck^ 

- [sKk'iz 






for IR region of kt 


dvr 


ci^/fkt 
■ In ( A f^zA'^ 


for large z 


Sy'CiCsTT VV Cl 
where A is the ultraviolet cutoff. 

Similarly, the following integral is also useful, 

sm{kzz)kz 


dkz 


TT 


akz + b + ckf 2cA 


(< 


-B-z _ -B+z 


)■ 


In general case where a = ci + C 2 k^, b = c^kf, for small kt, we obtain 


TT 


2cA 


—B—Z 


) - — 
Cl 


TT / 


e V <=1 




With this result, we can evaluate the following integral. 


27r 

(2vr)= 


/■A 

dkz / dktk 


sm{kzz)kz 


t^t ~ 


dktkt TT 
L (27r)2ci 


akz + b + ckt 
-J^kh 


e v-i 


for IR region of kt 


-- KtC V =1 tor large z 

dvrci 

1 1 


Svr^CiCs z 


(D4) 


(D5) 


(D6) 


(D7) 


1. long-range correlations of diagonal components 

In evaluating the long-range correlation functions of diagonal components fl29H . the fol¬ 
lowing expectation values with jSf" = (3i{z) ± f3i{0) are useful: 
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(cos(gz + /3,ti,2,3)) 
((5(a;) cos(g^ + /33=^)) 
(5(0)cos(gz + /3^)) 
((5(a;)(5(0)cos(gz + /3^)) 


cos qze~^^^ (D8) 

=Fsing2:((5(a;)/93(0))e“^^3±)/2 (J39^ 

— sing^((5(0)/93(a;))e“^^3±^/^, (DIO) 

cosgze-<^3D/2 [(5(a:)5(0)) T (5(a;)/33(0))(<5(0)/33(a;))], 

(Dll) 


where 


A{6{x)/3sm = A} dke^'^^6{k)/3s{-k)) = 


1 r d^k 

2] (^ 


e^^-^Ss3iO;k) 


2j (27r)3 

a^,2qA^T 1 

_ 

STiaQ iM'^ qz' 

{SiO)Psix)) = -{6{x)^sm: 


d^k . . 2iqA‘^ae.2kz T 

t sin [kzZj --^ tor x \\ z 


M2 a ;2 




±e'‘-‘)mm(k)) = 


d^k 


(D12) 

(D13) 


2A2 J (27r)3 


(1 ± cos kzz) 5*3(0; k) 


T 


d^k 1 ± cos kzZ 


2A2 J (27r)3 a;2 

oo 


In 


zA^ ’ 


(D14) 


167ra6,i 2^ 


(Atl2>/2 = '^<ik [{ft(*:)/?-(A^))(l ± e“-) ± ({I3,(k)l3;{k + 2qz))e‘’‘“ + cx.)] 

(1 ± COS kzz) 5 * 11 ( 0 ; k) ± ("sin A ;2 2:Ini5*i2(0; k) + cos A;^;2Re5*i2(0; 


1 1 

f d^k 

(1 ± cos k, 

- ± cos kzZ 

2A^J 
^ 1 

(27r)3 

f d^k ] 

2A^J 

✓ 

CO 

1 


1 00 




T 

_ Ivi AA 

^ IQiraQ^iA^u 

XXX f-k 

2q 


(D15) 
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From the above results, for instance, the {1,1} component in the direction reads 


fu{zz) = {{A+ S{z)){A + 6{0)) cos{qz + Psi^z)) cos 

X cos /32{z) cos/92(0) cos /3i{z) cos/9i(0)) 

= ^((^ + <^(^))(^ + <^(0)) [cos(g^ + /33^) + cos(gx + /33")]) 

x((cos/9^ + cos/5^) (cos/S}" + cos/9f)) 

~ [A^ cosgx (l — ((5 (x)/ 93(0))^) + 2Asingx((5(x)/93(0))l 

8 

~ -A^ cos (D16) 

8 

where we have used the fact that are logarithmically divergent, and terms including 

((5(x)/53(0)) oc z~^ drop faster than others for a large distance in the x direction. Also, 
((5(x)(5(0)) corresponds to a massive mode, which does not contribute to the long-range 
correlations. 

Similarly, in the x-y directions, we can obtain the following results: 


{5{xt) cosf3f) 

Pt) 


(5(0) cos 
(5(a;t)5(0) 


e 

0 , 

0 , 


-(/3±'>/2 


cos/3f) = e <^3±)/2 =F (5(a;t)/33(0))(5(0)/33(a;i))] 


(D17) 

(D18) 

(D19) 

(D20) 


where /Sf = l3i{xt) ± A(0), and 


(Ati,2,3)/2 = 


oo 




(D21) 
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